CURVES OF STEEPEST DESCENT ARE ENTROPY SOLUTIONS FOR A 
CLASS OF DEGENERATE CONVECTION-DIFFUSION EQUATIONS 
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Abstract. We consider a nonlinear degenerate convection-diffusion equation with inhomoge- 
neous convection and prove that its entropy solutions in the sense of Kruzkov are obtained as 
the — a posteriori unique — limit points of the JKO variational approximation scheme for an 
associated gradient flow in the L 2 -Wasserstein space. The equation lacks the necessary convex- 
ity properties which would allow to deduce wcll-poscdnoss of the initial value problem by the 
abstract theory of metric gradient flows. Instead, we prove the entropy inequality directly by 
variational methods and conclude uniqueness by doubling of the variables. 



1. Introduction 

The goal of this paper is to show on the example of the degenerate parabolic equation 

d t u = (u m ) yy + (b(y)u m ) y , y G M, t>0, (1) 

how the solution concepts of metric gradient flows and entropy solutions can be combined to 
obtain global in time well-posedness of the initial value problem. Specifically, we define an energy 
functional J 7 , construct time-discrete curves of steepest descent in the landscape of T with respect 
to the L 2 -Wasserstein metric by means of the JKO variational scheme [16|, and show that these 
curves converge in the limit of continuous time to entropy solutions a la Kruzkov [HI [20] (or, more 
precisely, in the sense of Carrillo [3J) for ([T]). This way, existence of solutions to (JXJ) is obtained by 
a variational method, and their uniqueness follows from the theory of entropy solutions. 

In {!]), the exponent m > 1 determines the degeneracy of the diffusion for vanishing densities, 
and 

b e L 1 (R) n W 1 '°°{M.) (2) 
is a given coefficient modelling heterogeneity in the convection part. As initial condition for (fT]), 
we prescribe a datum u° G i 1 (R) n L°°(R) of finite second moment. 

The link to metric gradient flows is established as follows. There is a spatial coordinate change 
y = T(x) and an according transformation p(t,x) = T'(x)u(t,T(x)) (see subsection 12.11 below) 
such that all sufficiently regular solutions u to ([1]) become solutions of 

d t p={p[a(x)p m - 1 ] x ) x , (3) 

and vice versa; here a £ W 2 '°°(M.) is a strictly positive function, determined from b via T. Well- 
known formal arguments, see e.g. ^7\, indicate that the time-dependent density functions p(t) 
satisfying ([3]) are "curves of steepest descent" in the energy landscape of the entropy functional 

F[p] :=- f a(x)p m dx 7 (4) 

with respect to the 2-Wasserstein distance. 

Indeed, if J- would be a geodesically A-convex functional, then the theory of A-contractive 
gradient flows — see e.g. [1] — could be applied to conclude the existence of a unique gradient 
flow for T in the space of probability measures. This flow's curves would be weak solutions to 
([3]), and by inversion of the coordinate transformation above, we could conclude well-posedness 
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for ([T]). Unfortunately, J- docs apparently not have the required convexity property; see Section [3] 
for details. 

Nevertheless, the variational structure behind (|3|) is on the basis for our proof of existence of 
solutions to (flj . Specifically, we employ the JKO (or "minimizing movement" ) approximation 
scheme to obtain time-discrete curves p T of steepest descent for J- ' . We prove that any family of 
such discrete approximations p T possesses a weak time-continuous limit curve p* . By inversion of 
the coordinate change, this provides a candidate u* for a solution to ([T]). We emphasize that this 
construction does not require geodesic A-convexity for T\ boundedness from below, coercivity and 
lower semi-continuity are sufficient. 

At this point, the variational framework of minimizing movements provides strong tools which 
allow us to prove that u* is actually an entropy solution for (flj in the sense of Carrillo jSj: we 
show that ti m € if oc (0, oo; 7? 1 (K)), and that for every non-negative test function ip € C^°(R + x R) 
and for every k 6 R + : 



We also obtain a non-trivial lower bound on the dissipation 2D(u„); see Proposition ^. 1 II for details. 

The derivation of the entropy inequality ([5]) is the core element of our proof. First, a time- 
discrete version of this estimate is proven directly for the JKO scheme, and this is passed to 
the limit. In the derivation of the discrete estimates, the key idea is — as usual — to choose 
appropriate variations of the minimizers. Here we build on the ideas that have already been 
employed in [16] for the derivation of the weak formulation of the linear Fokker-Planck equation 
and have been generalized later in [22] to the "flow interchange lemma" , see Lemma l4~2l variations 
are performed by an auxiliary gradient flow, which — in contrast to the gradient flow of ([J) itself 
— is A-convex and thus satisfies certain variational inequalities. 

Once that ([5]) has been establish, we adapt the doubling of the variables method in [18] to our 
case and show that u* is actually the unique entropy solution for the given initial condition u . A 
posteriori, we conclude uniqueness of the limit curve of the approximation scheme as well. 

The following theorem summarizes our main result in an informal way; the precise statement 
is given in Theorem 12.71 in Section 12.41 For a possible slight generalization avoiding the scaling, 
see Remark 12.91 

Theorem 1.1. Let an initial condition u° G L°°(K) of finite second moment be given. Then every 
curve obtained from the JKO approximation for T in the limit of continuous time corresponds - 
by a chance of coordinates — to the unique entropy solution for . 

A related observation about the connection between entropy solutions and gradient flows has 
been made recently by Gigli and Otto [15| in the context of the inviscid Burger's equation, see also 
previous results in [5] [21] ■ The interpretation of the coincidence between entropy solutions and 
gradient flows is that both types of solutions can be characterized by diminishing an underlying 
entropy functional "as fast as possible" . For previous results on the well-posedness of scalar 
conservation laws in Wasserstein spaces, we refer to [2] [4]. 

Apart from revealing an interesting connection between the two seemingly unrelated theories of 
entropy solutions and Wasserstein gradient flows, our simple example indicates a possible general 
strategy to prove existence and uniqueness of certain nonlinear evolution equations which can be 
cast in the form of a Wasserstein gradient flow of a functional which is not necessarily geodesically 
A-convex. First, use methods from the calculus of variations to construct a candidate for a solution; 
a priori, there might be several. Second, show that this candidate is an entropy solution by 
deriving further a priori estimates in the variational framework. Third, conclude uniqueness of 
the entropy solution. This strategy provides a new method to obtain entropy solutions, alternative 
e.g. to the classical vanishing viscosity approach [8], to the wave- front-tracking algorithm [7], or 
to semigroup theory [5]. In particular, the variational approach does not require the solution of 
auxiliary regularized problems. 




I |u* - k\<ft dydt 



(5) 




I Sgn(u„ - k)([{u™) y + b(u™ ~ k m )]tp y - fe y F'v) dydt > S)(tt,)b] > 
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To conclude this introduction, wc stress again that the goal of this paper is to establish a 
link between gradient flows and entropy solutions but not to refine the results on existence and 
uniqueness of weak solutions to ([1]) . This said, we remark that — as far as our model equation ([I} 
is concerned — uniqueness of weak ^-solutions could be proven also by other means, e.g. with 
the methods developed in [23] . 

The paper is organized as follows. In Section [2] we formulate the problem and state our 
main results. Section [3] recalls some basic fact from the theory of gradient flows. In Lemma 
13.41 we provide a result on the contractivity of quite general reaction-diffusion equations in the 
Wasserstein metric which might be of independent interest. In Section [4] we prove convergence of 
the scheme and show that the limit curves are entropy solutions. Section [S] contains (for the sake 
of completeness) the uniqueness proof of entropy solutions. 



Notation 

Measures, densities and Wasserstein distance. Here "p2(R) denotes the set of probability 
densities p € L 1 (R) with finite second moment / x 2 p(x) dx. Note that the symbol T-^R) is 
frequently used in the literature for the (wider) space of probability measures on R, which we 
shall denote by T-^R) instead. Given a measure p G 7^2 (R), its probability distribution function 
U : R — >• [0, 1] and pseudo-inverse G : [0, 1] — > R U {±oo}, respectively, are given by 

U(x) = fx((-oo,x)), G(w) = sup {a; G R | U(x) < u}. 
The space ^(R) is endowed with the L 2 -Wasserstein distance W2, defined by 

W 2 (ji,ji)=(J [GH-GH] 2 d w j 1 , (6) 

where G and G are the pseudo-inverse distribution functions of p and //, respectively. The pair 
(7^2 (R), W2) is a complete metric space. We refer to [27] for a more detailed explanation. 

Mollifications. We will frequently use the following elementary functions: the absolute value 
Abs(x) := |x|, the positive part PosPar(x) = (x)+ = (x + \x\)/2, the sign Sgn(x) = x/|x| (with 
Sgn(0) = 0), and the unit step function Stp(x) = (x + |x|)/2|x| (with Stp(0) = 0). In fact, we 
will mostly use their regularized versions obtained by mollification: denote by 5i : R — > R the 
standard mollificr 

Si{y) = |z- 1 cxp[-l/(l-y 2 )] for all ye (-1,1), 
1 otherwise, 

where Z > is chosen s.t. Si has unit integral. For e > 0, define the e-mollifier S e : R — > R by 
S e (y) = e~ 1 5i(e~ 1 y). Accordingly, we denote by Abs e = Abs*<5 e and PosPar e = PosPar*£ e the 
e- mollifications of Abs and PosPar, respectively. Notice that Sgn^ = 2<5 e and Stp^ = 6 e , which 
means in particular that Sgn e and Stp e are non-decreasing functions. 

2. Statement of the problem and results 

2.1. Coordinate transformation. In this subsection we shall establish the correspondence be- 
tween equations (TTJ and More precisely, wc prove the following proposition. 

Proposition 2.1. Let b 6 W 1 '°°(M.) n i x (R). Then, there exist a function a with the properties 

(al) a(x) > a > for all x G R, 
(ad) a G W 2, °°(M.), 

a bijective change of coordinates y = T(x) on R and a corresponding transformation S : L (R) — > 
i x (R) with 

S[u]{x) := T'(x)u(T(x)) (8) 

such that the transformation p(t) := S[u(t)] of an arbitrary weak solution u to ^ with initial 
datum uq G L 1 (R) ; satisfies equation ^ with initial datum p° = S[u ]. 
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Proof. Choose a positive constant a , and define 



a(y) = a exp ( - (m ^ f b( V )dr) ) ■ (9) 



2m 

The assumptions on b ensure a £ W 2 '°°(M.) and a > a for some a > 0. Therefore, the initial value 
problem 

% x (x) = a(1(x)) for x £ R, 1(0) = (10) 
admits a unique global solution % : R — > R, with X x (x) > a. By the chain rule, T satisfies further 

I xx (a;) = a / (X(i))T a ,(x), (11) 

and in view of we also have 

777 — 1 

a'o1 = -^— (ftol)^. (12) 
Now, let w be a weak solution to (JTJ), that is 

/•GO P pQO P pOO P 

- uip t dydt= / u m (p vv dydt- / / u m b<p y dydt (13) 

JO JO JR JO Jr 

for all test functions tp £ C£°(R+ x R). Let p(<) = S[«(t)] for all t > 0, which means that 

for all Further, for given <p, define the transformed test function ip £ C£°(R+ x R) by 

%l)(x,t) = <p(1(x),t). 

Using (TIT))) and (|11[) . one easily verifies that 

^(^(a;), t)X x (a;) = ip x (t,x), (p vv ( c I(x),t) c I x (x) 2 = ip xx (x,t) - ip x (x,t)a'(1(x)). (15) 

Substitute (fT4| into (fT3|) . perform the change of variables (y, t) = (1(x),t) under each of the 
integrals, and simplify the expressions containing test functions using (|15l) . This yields 

/oo p 
I p m ^ xx %- {m+1) dxdt 

(16) 



Defining a : R -> R by 



^ y p m ^[(a'oT)T;( m+1 ' + (6o!r)^ m ] dedf. 



a(as) = 7 (aoI(i))-(™ +1 ', (17) 

m — 1 



we find by direct calculations that 

l ax W -H^±±(a' o T)T-(— « # (a' o I)!^"* 1 ) + (6 o T)T~ m . (18) 

TO TO — 1 

Substitute (JTTJ) and ([T8|). respectively, in the first and the second integral on the right hand side 
of (US), to find 



oo f \ r°° r l 

p^t dx dt = / / ap m ijj xx dxdt / / a x p m ip x dx dt. 



i) 



It is easily checked that this is a weak formulation of ([3]). □ 



Remark 2.2. The function a in Provosition \2. 1\ is not uniquely determined: different choices of 
a > in §§§ change a by a positive factor. On the other hand, if the function a is given, then 
b and the corresponding change of variable y = T(x) can be recovered from a in a unique way. 
Indeed, it follows from (jlOp and (|17p that 



1(x) 



f 


~m - 1 ' 


10 


TO 



'■■I 1 



d£, (19) 



ENTROPY SOLUTIONS VIA WASSERSTEIN GRADIENT FLOW 



5 



while (|17p and (|12p imply 



2m 

bot=^—(loga) x . (20) 



Remark 2.3. The scaling X is still well-defined if the assumption b G L (R) is relaxed to b oc |y| _1 
as \y\ — > +oo. However, assumptions (al) & (a2) may be not satisfied in those cases. Therefore 
we require b G L 1 (M). Let us also emphasise that the case b = const cannot be included, since the 
solution to would blow up at some negative x. 

2.2. Entropy solutions for The notion of entropy solution for ([l) used here is a variant of 
the one originally introduced by Carrillo [3] and later adapted by Karlsen et al. [T5]. To motivate 
this definition, consider the usual viscous approximation of ([1}, 

d t U = (b(y)u m )y + (u m ) yy +uu y y, (21) 

which possesses smooth and classical solutions u u for every v > 0, provided that the initial 
condition u° is smooth enough. The following formal considerations arc made under the hypothesis 
that it„ converges — locally uniformly on R + x R and in Lf oc ([0, T[xR) — to a limit function as 

Similarly as in the classical approach by Kruzkov |20j for scalar conservation laws, we would 
like to derive from (|21j) an evolution inequality for all functions of the form \u — k\ with given 
k G R+. As usual, the calculations are carried out with a suitable approximation of \u — k\; in 
our case, we multiply equation (f2"Tj) by Sgn e (u m — k m ) with a mollification parameter e > and 
rewrite it in the following way: 

Sgn £ (u m - k m )d t u = Sgn £ (u m - k m ) [{u m ) y + b(u m - k rn ) + vu y \ y + Sgn e (u m - k m )b y k m . (22) 

Now let T > 0, integrate ([22]) against a non- negative test function (p G C^°(]0, T[xR), and integrate 
by parts in the first term on the right-hand side: 



if Sgn e (u Tn - k m )d t u dy dt = (23) 

Sgn> m - k m )([{u m ) y + b(u m - k m ) + uu y ]tpy - b y k m ^ dydt (24) 

j <pbSgn' e (u m - k m )(u m - k m )(u m ) y dydt (25) 

+ f f Sg^ e {u m -k m )[{u m )l+mvu m - l ul]^dydt. (26) 
For further simplification, observe that inside the integrand in (|23j). 



limSgn e (w m - k m )d t u = Sgn(u - k)d t u = d t \u - k\ 

at every point (i, y) with u(t, y) 7^ k. Likewise, in (|24[) . the term with v becomes Sgn(w — k)u y = 
\u — k\ y in the limit e \. 0. After integration by parts, its contribution is 

-v i \u- k\ifyy dydt. 
J a Jr 

This term is negligible in the limit v \. 0. Next, the integrand in (|25|) can be rewritten as 
(pb[R € (u m - k m )] y with R e (s)= f rSgn' e (r)dr. 



Since R e converges to zero uniformly for e J, 0, it follows — after an integration by parts — that 
the integral (|25|) vanishes in that limit. Finally, the integral ([26]) obviously gives a non-negative 
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contribution. Thus in the limit e 1 0, and after an integration by parts with respect to time in 
(|2"3"|) , the equality ([2"3" ]) -([2l) ]l implies that 

f [ \u-k\<p t dydt> [ [ Sgn(u-k)([(u m ) y + b(u m -k m )]<p y -b y k m <p)dydt 
Jo Jr Jo Jr v ' l07 \ 

+ limsup / / Sgi<(u m -fc m )[(u m ) y ]Vdydt. 
We take this as the defining inequality for entropy solutions. 

Definition 2.4 (Definition of entropy solution). Let u° g L 1 ni°°(K).v4 non-negative measurable 
function u : R + x 1 1 is on entropy solution to ([T]) urai/i initial condition u° if u £ L 1 PI 
L°°(]0,T[xR) for all T > 0, if u m G £? oc (°> oo; ff x (M)), -> w° in L 1 (R) as f J, 0, and i/ 

inequality (|27|l is satisfied for all nonnegative test functions ip G C ( ?°(R+ x R), and for all k G R + . 

Remark 2.5. T/ie inclusion of the very particular dissipation term in inequality (|2T[) seems a bit 
ad hoc. This term, however, plays a key role in the proof of uniqueness, see Section^ Moreover, 
the choice of this term is less arbitrary than it appears: any smooth and uniformly convergent 
approximation of the sign function could be used in place of the mollification Sgn e there. 

Finally, we remark that a substantial part of the article |18j is devoted to proving that all 
functions u of a certain regularity which satisfy (|27[) without the dissipation term actually satisfy 
it also with the dissipation term. For us, the dissipation term results very naturally from our 
construction of solutions u. 

Remark 2.6. For a function u of the specified regularity, in order to be an entropy solution it is 
sufficient that inequality (|27[) is satisfied for all test function tp of the form tp(t,y) = 6(t)<j>(y) with 
arbitrary non-negative functions 9 G C^°(R+) and <j) 6 C£°(R). This follows immediately since 
the latter products lie dense in C^°(R+ x R), see e.g. [141 Theorem 4.3.1]. 

In difference to [18j . we do not require — and, in fact, cannot prove — that u is a continuous 
curve in L 1 (R) for t > 0. 

2.3. The JKO scheme. The JKO scheme [16] is a variant of the time-discrete implicit Euler 
approximation for the solution of gradient flows in the non-smooth metric setting of the L 2 - 
Wasserstein distance. In its core, it is a special case of De Giorgi's minimizing movement scheme 
[10) ; see the book [T] for an extensive theory. 

We apply the JKO scheme to the functional T defined in ((4]) with a defined implicitly in 
Proposition 12.11 and inital condition p° := S[u°], sec (JSJ. To this end, let a time step t > be 
given. For every a € ^(R), introduce the associated Yoshida penalization a) of T by 

Let further an initial condition p° € 7^2 (R) with T(p°) < +oo be given. We define a sequence of 
densities S 7-2 (R) inductively as follows: 

(1) p°:=p°. 

(2) For n > 1, let e T^R) be the (unique global) minimizer of T T {-; p" 1 ). 

In Lemma 14.11 we prove wcll-dcfinitcncss of this scheme. In the following, we denote by p T : 
[0,oo[— > ^(R) the piecewise constant interpolation of the sequence (p")„ 6 N, with 

p T (t) = p n T for (n - 1)t < t < nr. (28) 

2.4. Main results. We are now in the position to give the precise statement of the two main 
results of this paper. 

Theorem 2.7. Let p° g 7 , 2(R)ni°°(R), and let a satisfy conditions (al)&(a"2) of Proposition^^ 
Define discrete curves p T : [0,oo[— > ^(R) by means of the JKO scheme from Section \2.3\ above. 
Then every vanishing sequence {Tk)k£N of time steps contains a subsequence (not relabeled) such 
that the p Tk converge — in L m (]0, T[xR), and also uniformly in W2 on each time interval [0,T] 
- to a curve p* : [0, 00 [— > T-^R) that is continuous with respect to W2. The rescaled function 
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U* — S is an entropy solution to (JTJ) in the sense of Definition \2~~$\ with initial condition 

«° = s-V]. 

The proof of Theorem 12.71 is performed in Section 2J 

Theorem 2.8. For each initial condition u° £ V2 (R) H L°°(M.), there exists only one entropy 
solution to (JlJ in the sense of Definition \2~4\ Consequently, every sequence of time- discrete ap- 
proximations p Tk converges to the same limit . 

The proof of Theorem 12.81 is taken from [18], with minor modifications. In order to make 
the present paper self-contained, we outline the argument and review the relevant calculations in 
Section [3J 

Remark 2.9. The notion of entropy solution can be posed as well for the equation ^ without 
passing through the scaling established in Proposition \2.1l and a uniqueness theorem in the spirit 
of Theorem \2.8\ can be derived with arguments in [17] . Nevertheless, we opted for developing our 
theory based on the notion of entropy solution for the equation which can be seen as a porous 
medium equation with a nonlinear convection perturbation, and it is therefore of interest in the 
applications. 



3. Convex functionals and contractive gradient flows 

Contractive gradient flows (or "k-Hows" , see Dcfinition l3.1l below) constitute our key tool for the 
derivation of a priori estimates on solutions p to ([3]), or, equivalcntly, for solutions u to ([IJ. First, 
we recall the definition of K-flows and their relation to gcodesically A-convex functionals. Then 
we prove K-contractivity of a class of flows that is relevant for our needs. We refer to [1] [27, 28 
as references on the general gradient flow theory on Wasscrstein spaces. 

3.1. Fundamental definitions and relations. A curve p : J — > T-^R) is called absolutely 
continuous in W2 on the interval / C R if there exists a function g £ L\ oc (I) such that 

W 2 (p(t), p{s)) < j g(r)dT for all t, s E I. 

An absolutely continuous curve p : [0, 1] — > 7^2 (R) is a constant speed geodesic if 
W 2 (p(s),p(t)) = \t-s\W 2 (p(Q),p(l)) for all f, 8 6 [0,1]. 

Definition 3.1 (re-flow). A semigroup 6* : [0, oo[x"P2(R) — > ^(R) is a n-How for afunctional 
^ : ^(R) — > R U {+00} with respect to W2 if, for arbitrary p G ^(R)? the curve s H- ©|,[/o] is 
absolutely continuous on [0, oof and satisfies the evolution variational inequality (EVI) 

w 2 (eup],p) 2 + ^2(e%[p},p) 2 <nv)-n&l[p}) (29) 



2 da 



2 



for all s > 0, with respect to every comparison measure p <= 7^2 (R) for which ^(p) < 00. 

Remark 3.2. The symbol d + /dcr stands for the limit superior of the respective difference quo- 
tients, and equals to the derivative if the latter exists. 

The fact that a functional admits a K-flow is equivalent to the A-convexity of \& along 
geodesies. The characterization of K-flows by convexity will not play a role in our further consid- 
erations, but we cite the respective result for the sake of completeness, see [T] for further details. 

Theorem 3.3. Assume that the functional : ^(R) — > RU {+00} is A-convex (along geodesies), 
with a modulus of convexity A £R. That is, along every constant speed geodesic p : [0, 1] — > T-^R), 

n P (t)} < (i-tMp(o)}+m P (i)}-~t(i-t)w 2 (p(o),p(i)) 2 (so) 

holds for every t G [0, 1]. Then possesses a uniquely determined K-flow, with some K > A. On 
the other hand, if a functional possesses a K-flow, then it is X-convex, with some A > k. 
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3.2. A special class of K-flows. In this section, wc derive a sufficient condition under which the 
functional 

*(?/) = / F(x,r](x))dx (31) 
Jr 

for a sufficiently smooth function F:Ix R + — ^ R admits a k-How. For a concise formulation of 
that condition, we introduce the adjoint function H : R x R + — ► R by 

£T(a:,0=^(a:,l/0. ( 32 ) 

which satisfies 

H 3t (x,Z)=tF x (x,l/t), Ht(x,Z)=F(x,l/t)-l/ZF v (x,l/t). (33) 

Lemma 3.4. Lei F e C 2 (l x R + ) be given, and assume that there exist constants < c < C 
such that 

c<r)F vv (x,n) <C and \F xr ,(x, T))\ < C for all (x , r/) g R x R+ . (34) 
.Assume further that there is some k S R smc/i i/iat 

(^O^iJ^O-^ 2 (35) 
is (jointly) convex on R x R + . Then, the solution operator of the evolution equation 

dtV = 'D x (r)D x [F v (x,ri)}) (36) 
is a K-flow for the functional \& /rom (|31|. 

The convexity condition (|35|) can be rephrased as 

-.2 rr _ f H xx - K H X £ 



D H = „ rr i s positive semi-definite for all £), (37) 



and the condition (|34|) is equivalent to 

c < £ 2 i%(a;,£) < C and |Lf x (z,£) - £-H"a*(a;,£)l < C for all a; g R, £ g R+. (38) 
Proof. We need to show that the EVI 

i^w 2 (©U ^) 2 + \ w 2 (eW,v) 2 < m - *(®W) (39) 

holds for all smooth and strictly positive densities 77 , 77 g V<x (R) at almost every time i > 0. 

Instead of proving (|39j) directly for we prove it for a family of regularizations and then 
pass to the limit. For every N £ N, the regularized functional vffjv is defined by ^^(ry) = ^(77) 
if 77 e ^(R) is supported in [— N, N], and ^n(v) = +°° otherwise. We show that the solution 
operator &n to the boundary value problem 

dtV = T) x (iqF m {x,ri)ri x +r}F vx (x,ri)), n x {t,N) = r) x (t,-N) = (40) 

is a K-flow for i.e., the analogue of (|39p holds. The densities 77 and rf are approximated by 

Vn = (v + Ajv)1[-jv,jv]j Vn = i 7 f + jv,jv] > 

with constants /jjv, 7^ > such that 77 jv, 77^ £ ^(K). Notice that ^jv(^jv), ^w^jy-) converge to 
^(77), ^(r; ), respectively, as JV ^ 00. 

By the lower bound on r\F m required in (|38|) . the equation (|40[) is uniformly parabolic. Thus the 
boundary value problem with initial condition 77^ possesses a solution 77 with 77(f) = &%{r)%) for all 
t > such that each 77(7;) g T> 2 restricts to a smooth and strictly positive function on [-N, N] and 
vanishes outside of that interval. The associated distribution function U : [0, 00 [x [— N, N] — > [0, 1] 
with 

nX 

U{t;x) = / v(t;y)dy 
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is a smooth diffeomorphism from [—AT, N] to [0, 1], for every t £ K+; this follows from smoothness, 
strict positivity and mass preservation of rj. The inverse G : [0,oo[x[0, 1] — > [— N, N] of U with 
respect to x satisfies by definition 

G(t;U(t;x))=x for all x £ \-N,N}. (41) 

In particular, we have 

G(t; 0) = —N and G(t; 1) = N. (42) 

We differentiate equation (|4ip with respect to x and with respect to t, respectively, to obtain 

l=r)(t;x)d u G(t;U(t;x)), and (43) 
= d t G(t; U(t; x)) + d t U{t; x) d u G{t; U(t; x)). (44) 

On the other hand, integration of (ptfj)) with respect to x yields 

d t U(t; x) = r](t; x)B x [F v (x, rj(t; x))] . 
We multiply by d u G(t; U(t; x)), substitute (gj) for d t U and (gSJ) for V- This yields 

d t G(t;U(t;x)) = -B x [F n (x,ri(t;x))}. (45) 

We rewrite this (omitting the dependence on t) using again ([4"3")) and the properties of H from 
(33): 

-D,[F„(x,77(x))] = —^—-(p x [F(x,r)(x)) — r](x)F v (x,r](x))] — D x F(x,r)(x)) + d x r](x)F v (x,r}(x))) 
= -}-T> x [Hs(x, 1/ V (x))} - H x (x, 1/ V (x)) 

= d u G(U(x))T) x [Hz(G(U(x)),Gu(U{x)))} - H x (G(U(x)),d u G(U(x))). 
In combination with (|45[) . we have 

d t G = T>„[Hz(G,d u G)} - H x {G,d^G). (46) 
Further, observe that a change of variables to = U(x) leads to 

<M*) = /_%(*, ^))dx = f%{G(U(x)), ^L^) dx 

f d u G{u)F(G(u), _ * , )du= f H(G(u),d u G(u))du. (47) 
o V d u G(u)J J 

Now, let G be the inverse distribution function of 77 jy. Using the representation ([6]) of the Wasser- 
stcin distance, we obtain 

W 2 (77(i),r?) 2 = / [G(t;w)-G( W )] 2 d W . (48) 
Jo 

Now we can prove (|39|) . Combining (j48|) with the evolution equation (j46|) and the representation 
(03 of we find 



~W2(ri(t),VN) 2 + ^W 2 (r 1 (t),f) N ) 2 = G t [G-G]dw + |^ [G-G] 2 d^ 

= / D u [H e {G,G u )][G-G]du- f H X (G, G U )[G - G] da; + £ / [G-G] 2 dw 
Jo Jo ^ Jo 

= } / ff ? (G,G a) )[G w -G w ]dw + / fl a (G,G w )[G-G]dw + ^ / [G-G] 2 d^ 
Jo Jo ^ Jo 

(**) f 1 ~ - 

< / [H(G,G u )-H(G,G u )]du = * N (fi N )-* N (ri(t)). 
Jo 



in 
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In step (**) we have used the joint convexity of (pj|) . The equality (*) is established integrating 
by parts; the boundary terms vanish because of the boundary conditions (|4"2l satisfied by G and 
G. This proves the analogue of (|39[) for each VP^r with its respective flow ©at. 

Finally, we pass to the limit N — > oo. For definiteness of notation, we denote the solutions to 
(|40|) by r]N(t) = & N (r] N ) from now on. 

To begin with, we show that a subsequence (ryjv) C (?7jv) converges to a limit function 7700 
satisfying the unregularized evolution equation (|3"6")) . For the following calculations, the bounds 
required in (|34[) arc important. First, we derive an ^(HQ-estimate: 



T)% dx = - rjNFrmiVN, x) (t?jv)^ da; - / {vn) x F xri (n N , x)rj N dx 

-N J-N J-N 



N r N r 2 

2 J„ C / r„_ A2 j_ , ° / „2 



< -cy O^)* dT - 2 y (^)x da; + y ^ dx. 

Integration with respect to t £]0,T[ yields an A^-uniform bound on t)n in L 2 (0,T; H 1 (M.)) in 
terms of 1 1 77 1 1 ^2 . By Alaoglu's theorem, this provides weak convergence of a subsequence t)n> to 
a limit r?oo m L 2 (0, 00; Second, we prove an i/ _1 (R)-estimate on dtf]N- To this end, let 

®xx = d t {rjN) w hh 9 x (t,N) = 9 x (t, —N) = 0, and calculate, using again (j3~4"|) : 

9 2 x dx = - / 6*6^ dx = - I 6d t (r} N ) dx = - 9 X (r] N ) x rj N F m dx + 6 X 7] N F VX dx 

-N J-N J-N J-N J-N 

1 rN pN pN 

<- eldx + c 2 ( m )l dx + c 2 V 2 N dx, 

L J-N J-N J-N 

which, in combination with the H 1 (R)-bound from above, shows that dtijN is uniformly bounded 
in L 2 (0, T; H _1 (R)). The Aubin-Lions lemma, see e.g. [26], shows that the convergence of 777V to 
7700 is actually strong in Lf oc (M. + x R). 

Strong convergence in Lf oc (M.+ x R) and weak convergence in L 2 (0, 00; i/ 1 (M)) allows us to pass 
to the limit in the weak formulation of P0|) . 

(p t r)N> dxdt = - / / (p x [T)N'Fr rn (x,r]N')(T)N')x + 'nN'Fxri(x,r]N')]^ xdt 



where cp £ C£°(R+ x R) is a test function. Here we use that the nonlincarities r)F m (x,ij) and 
r)F xri (x,r]) grow at most linearly with respect to 77 by assumption (|34[) , and so they converge in 
L 2 3C (R+ x R), too, by dominated convergence. We conclude that the limit 7700 is a weak (and thus 
— by parabolic regularity theory — also the unique classical) solution to the Cauchy problem (|36|) 
with initial condition 77 . In other words, we have 6^,(770) = 7700(7;) for almost every t > 0, where 
(Sip is the solution operator to (|3"6"|) . 

Passage to the limit in the EVI ([39]) is more subtle because the convergence of the time deriva- 
tives dtrjN' is a priori only weak in £ 2 (0, 00; iJ _1 (R)). In [9], an equivalent time-integrated version 
of (f39| is given, and the latter is stable under strong L 2 oc - convergence. We refer the interested 
reader to [9J Section 3] for further details. □ 

Lemma 13.41 gives an indication that the functional J- from (jj) is not geodcsically A-convex. 
Note that, by Theorem 13.31 convexity of J- would imply the existence of a /c-flow for fy. However, 
with F(x,r)) = a(x)rj m , we have H(x,^) = a(a;)£ 1-Tn , and thus 



D 2 H(x,t) = 



a"(x)^- m - k -(m - l)a'(x)£- m 
-(m - l)a'(x)£- m m(m - l)a(ar)^( m+1 ) 



Unless a is a constant, there must exist an x £ R with a"(x) < 0. Thus, no matter how n is 
chosen, there exists further a sufficiently small £ £ R+ such that the top left element of D 2 H(x, £) 
is negative, that is a"(x)£ _ (" i_1 ^ — k < 0. Consequently, D 2 H(x,£) is not positive semi-definite. 
Note, however, that we have only shown sufficiency of the joint convexity of (|35j) for the existence 
of a K-flow. 
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4. Construction of entropy solutions 
This section is devoted to the proof of Theorem 12.71 

Throughout this section, p™ denotes the nth iterate in the JKO scheme for given time step 
r > as described in Section [2~3l see also (|49|) below. Further, p T : [0, oo[— > Vs(M.) denotes the 
piecewise constant interpolation of the sequence (p™) n6 N, with (|28[) . 

4.1. Well posedness of the JKO scheme. Recall the inductive definition of the p™, 

p n T £ P 2 (R) is the unique minimizcr of 7V(p; p™~ x ) = -^W 2 (p, p™" 1 ) 2 + F{p). (49) 

2r 

Lemma 4.1. The time- discrete scheme is well-defined in the sense that for any given p" -1 £ 
V%(^), the functional J>(s Pt ) indeed admits a unique minimizer p™ £ ^(R). 

Proof. We argue that the direct methods from the calculus of variations apply to the minimization 
problem for T r (•; p?" 1 ) on V 2 (R). 

First, observe that p" _1 £ T-^R) implies that the J r (p"^ 1 )-sublcvel of this functional is non- 
empty and weakly-* relatively compact. Indeed, p"~ x itself belongs to that sublevel, and if p* £ 
T-^R) is any other element of this sublevel, then W 2 (p*,p™ -1 ) 2 < J^p™ -1 ), since J 7 is a non- 
negative functional. This bound on the Wasserstein distance provides a control on the second 
moment of p* in terms of p™ _1 alone, and so the sublevel is a tight collection of measures, and 
thus relatively compact in the weak-* topology. Further, both T and the distance W 2 (-, p™ -1 ) 
are lower semi-continuous with respect to narrow convergence. This, together with the fact that 
J r T (-; p™ 1 ) is bounded from below, guarantees the existence of at least one minimizer p™ £ ^(R). 
Finally, since necessarily -^(p™) is finite, p" belongs to the subspace of densities "P 2 (R). 

To prove uniqueness of the minimizer, it suffices to observe that the squared distance W 2 (-, p" -1 ) 2 
is convex, and that the potential JF is strictly convex (recall that m > 1) in the sense of linear 
interpolation of measures. Thus, there can be at most one minimizer. □ 

4.2. Basic energy estimates and flow interchange. Immediately from the construction in 
(|49|) . we obtain the canonical energy estimate 

1 N 

— W 2 (p?, pr 1 ? < HPr) - HPr) for all N £ N. 

n— 1 

which induces the following uniform bound on in L m (]R): 

Hp?\\™ m <F(Pr)<F(P°)- (50) 

m 

Since T is non-negative, it follows further that 

oo 

^W 2 (p;\pr 1 ) 2 <2r^(p?). (51) 

71=1 

For our purposes, we derive stronger a priori estimates by variations of the minimizers p™ along 
specific k- flows, following the general strategy from [22l Section 3]. 

Lemma 4.2. Let \I r : ^(R) — >] — oo, +oo] be a lower semi- continuous functional on T-^R) which 
possesses a n-flow Define further the dissipation of T along by 

D*(p) := limsup-LF(p) - J"(6^p)l 

for every p £ ^(R). // p™ a«c? p™ are iwo consecutive steps of the minimizing movement scheme 
flM]), then 

v^p"- 1 ) - *(p») > rS)*(p?) + |w 2 (p»,pr 1 ) 2 . (52) 
In particular, ^(p"" 1 ) < oo implies 2)^(p") < oo. 
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Proof. Since (|52p is trivial for vE^p" -1 ) = +00, there is no loss in generality to assume ^(p™ -1 ) < 
00. We can thus use the EVI inequality (|29[) with p := p™ for the comparison measure p := p" -1 . 
By lower semi-continuity of we conclude 

^(pr 1 ) - > limsup (^(pr 1 ) - *(6^p; 1 )) (53) 

sj.0 

w 2 (©^,pr 1 ) 2 ) + Jw 2 (p^,pr 1 ) 2 . (54) 



1,. / d+ 
> — lim sup I 

2 sj.o 



V da 



lim sup 

sj.0 



W 2 (6^p^,pr 1 ) 2 ) >limsupi(W 2 (6^p",pr 1 ) 2 -W 2 (p^pr 1 ) 2 ). 



Notice that we have used the (absolute) W 2 - continuity of the curve s H> 6|,p™ at s = in the 
first step. The absolute continuity implies further that 

d+ 
da 

By definition of p™ as a minimizer of J- T from (|49|) , we have 

for every s > 0, and thus 

W 2 (6^?,rfr 1 ) 2 - W 2 (p>r') 2 > 2t[F{^) - H<S%Pr)}- (56) 
Inserting flSgJl into (|55]). and then (J55J) into (JM) yields ([52]). □ 



Corollary 4.3. Under the hypotheses of Lemma \4-2\ let the K-flow 6$ be such that for every 
n£N, £/ie curve s 1— > ©^p" Zies in L" l (R), where it is differentiable for s > and continuous at 
s = 0. And moreover, let a functional M. : P 2 (M) — ►] — 00, 00] satisfy 



liminf ( — 



T{<5% P n T ) >«)■ (57) 



s|0 V dcr 
T7ien </ie following two estimates hold. 

For every n € N: ^(p?- 1 ) - > T.£(p; 1 ) + |\V 2 (p™, p?" !) 2 ; (58) 

N 

for every N e N: *(pf ) < *(p°) - r ^ £(p") + r max(0, -k)F( P °). (59) 

Ti=l 

Proof. First, wc estimate 2)^ from below by .ft. The L m (R)-regularity assumptions on 6* imply 
that s h-> J-"(6*p") is differentiable for s > and continuous at s = 0. Thus, by the fundamental 
theorem of calculus, 



D*(p; 1 ) = hmsup i[^(p») - F{&lPr)] = limsup f ( - -f -F(©*p; 1 )) dz 





d 


/ lim inf ( — 




Jo s io V 


~da~ 



by Fatou's lemma and ([57]) . Now estimate ([58)1 follows directly from ([52]) . Estimate (|59|) is obtained 
by arranging (|58[) in a telescopic sum and combining it with the energy estimate (|51[) . □ 

4.3. Refined a priori estimates. The main ingredient in order to get suitable compactness is 
the following a priori estimate. 

Lemma 4.4. There is a constant A depending only on p° (and in particular not on t) such that 
the piecewise constant interpolants p T satisfy 

\\P? /2 \\mo,T;m m )- A{1 + T) f° rallT>0 - ( 6 °) 
In particular, p T (t) m/2 £ ff^R) for every t > 0. 

For the proof of Lemma 14.4) we need a r- uniform control on the second moment of p r , which 
is established in a further lemma. 
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Lemma 4.5. There is a constant B depending only on p° such that 

[ \x\ 2 p T (T,x)dx < B(l +T) forallT>0. (61) 



The proof of Lemma 14751 is obtained by a straight- forward application of the "flow interchange" 
Lemma 14.21 Specifically, we choose as auxiliary functional ^ := m the second moment, 



m (p) = / x p( x ) 

Jr 

This functional possesses a K-flow (with k = 2), which is explicitly given by 

(e s mV )(x)=e 2s V (e 2s x). 

Clearly, s M- &^r] defines a curve in L m (R) for every r\ G L m (R); this curve is differentiable for 
every s > 0. 



Proof of Lemma \4-5\ In order to obtain a sensible bound M. m in (|57[) , we compute the s-derivative 
of T along 6 s m for a given 770 e 7> 2 (R) H L m (M): 



d 



s=o v ds 



s=0 



y ds s=o \ Jr m ) 



2(m-l) / — 770(2) dz- 2 / -i;o(zdz<MJ(tjo) 



with the constant 

A/ := 2(m - 1) - 2 inf ^$ > 0. 

zeR a(z) 

We apply Corollary 14.31 with R m := —A/ J 7 ; estimate (|59|) becomes 

/ x 2 p"(x)dx< f x 2 p°(^)dx + r V £ m (Pr) < [ x 2 p°(x) dx + MF{p°) Nt, 
Jr Jr n _ 1 Jr 

using the monotonicity ([50]) of the potential. From here, (jBTj) follows immediately. □ 



We now turn to the proof of Lemma 14.41 above. Also here, the key element is an application of 
Lemma 14.21 This time, the auxiliary functional := H is the entropy, 



H (p) = / p{.x)\ogp{x)dx. 
Jr 



One of the celebrated results from the theory of optimal transportation is that this functional 
possesses a /s-flow 6h, with k = 0, which is given by the heat semi-group, see e.g. [§1 EH [2T| . 
More precisely: for given 770 <G the curve s H> rj(s) := ©^770 s °l ves the initial value problem 

d s ri = i] X x, 7?(0) =770, (62) 

in the classical sense: 77(5) is a positive density function for every s > 0, it is continuously 
differentiable as a map from R + to C°°(K) n L 1 (M), and if 770 <G L rn (R), then r](s) converges to 770 
in L rn (R) as s i 0. 

In order to extract information from the flow interchange estimate for H, we need another 
technical ingredient. 

Lemma 4.6. There is a constant C depending only on p° such that 

-C(l + T)< H(p T (t)) < C for all T > 0. (63) 

Proof of Lemma \4-6\ One verifies by elementary calculations that 

2 1 

— s 1 / 2 < slogs < — s m for all s > 0. 

e (m — l)e 
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And so, for every r\ € ^(ffiQ, we have on one hand that 

Hfa) = f r,(x)logr,(x)dx < 1 / »,(*)"* da; < m 

Jk (m-l)eJ R a(m-l)e 

while on the other hand, 

H(r?)>-- / (l + ^r^l + x 2 ) 1 / 2 ^)- 1 /^ 



>--(/ T^V^f /( 1 + ^)^)^V /2 >-^f 1 + f ^r}{x)dx 



e V ./» 1 + x- 



Now, with p T {t) in place of rj, it is straightforward to conclude both estimates in by using the 
energy bound (|50p and the moment control (joTj) . respectively. □ 



Proof of Lemma \4-4\ By the preceding discussion, ©h satisfies the hypotheses of Corollary | 
we shall define a suitable lower bound .%i for the use in (|57|) . By the spatial regularity of r](s) for 
every s > 0, the following calculations are justified: 

3/(6^0) = — I a(x)d s (n(s,x) m ) dx = I a(x)r)(s, x) m ~ 1 n xx (s, x) dx 



-(m — 1) / a(x)n(s,x) m 2 r] x (s,x) 2 dx - / a x (x)r](s, x) m 1 r) x (s,x)dx 
Jr Jr 

- 4(m ~ 1) /" a(x)[a x (r ? ( S ,x) m / 2 )] 2 dx+- / ^(^(s,!)"^ 



<- li! ^ / [^(^,xr /2 )] 2 d, + ^ f V ( S ,xrdx, 

771 Jr m Jr 

where we have used the notation := sup a xx < +00 in view of the property (a2) in Proposition 
12.11 Consequently, we define 

■MP) := 4(m ~ 2 1} - / [^(P m/2 )] 2 d^- ^ / P m dx. (65) 
TO Jm m Jr 

It is easily seen that (|57|) is satisfied with this choice of .ft: indeed, since ©fj is continuous in 
L m (R) at s = 0, it suffices to observe that Ku is lower semi-continuous in L m (M.) by Lemma [A. 11 
So Corollary 14.31 is applicable, and from (|5^|) . we conclude that 

n— 1 n— 1 

for every iV <G N. In combination with the bound ([S3"| on H, it follows that 

., JV JV lX N 

} --E ll(^) m/a &. * H ^°) " H (^) + — Iklir- + fc^r ^ 



r IlL" 

n— 1 n— 1 n— 1 



* — ; m 

4(m - 1) 



< C(2 + TVr) + 



T( P °)Nt. 



The energy estimate (|50|) has been used to derive the last line. From here, it is immediate to 
conclude (I5D1). □ 



4.4. Weak and strong convergence of the discrete curves. Our first convergence result is 
a standard consequence of the energy estimate. 

Lemma 4.7. Every vanishing sequence (rk)keN of time steps tu > contains a (non-relabelled) 
subsequence such that p Tk converges — uniformly on compact time intervals — in W2 to a Holder 
continuous limit curve p„ : [0, 00 [— > ^(K)- 
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Proof. From (f5Tj) . we conclude the Holder- type estimate 



W 2 (p T (t),p T (s)) < v^^ymax^li-sl) 1 / 2 (67) 

for arbitrary s,t > 0. The claim is now obtained as a consequence of (a refined version of) the 
Ascoli-Arzela theorem p] Proposition 3.3.1]. □ 

In the following, let a sequence (r^) be fixed along which p Tk converges to some limit p» in the 
sense described in Lemma l4~7l For the derivation of the entropy formulation (|27|) for it* = S _1 [ ( o»] 
in the next section, the following stronger convergence result is needed. 

Proposition 4.8 (Strong i m -compactness). The curves p Tk : [0,oo[— > T-^R) converge to p* in 
L m (]0,T[xR) for every T > 0. 

The proof of Proposition 14.81 is obtained via an extension of the Aubin-Lions lemma as given 
in [25]. The precise statement is recalled here for convenience: 

Theorem 4.9 (Theorem 2 in [IS]). On a Banach space X, let be given 

• a normal coercive integrand 3 : X — > [0, +oo], i.e., $ is lower semi- continuous and its 
sublevels are relatively compact in X ; 

• a pseudo-distance j:IxI-> [0,+oo], i.e., g is lower semi- continuous, and g{p,if) = 
for any p,r\ G X with $(p) < oo, $(r~i) < oo implies p = n. 

Let further U be a set of measurable functions u :]0,T[— > X, with a fixed T > 0. Under the 
hypotheses that 

pT r-T-h 

SU P / $( u {t)) di < oo and lim sup / g(u(t + h),u(t)) dt = 0, (68) 

U contains an infinite sequence {u n ) ne jq that converges in measure (with respect to t €]0,T[) to a 
limit it* :]0,T[— > X. 



Proof of Proposition \4-8\ Fix some T > 0. We verify the hypotheses of Theorem 14.91 for a specific 
choice of X, $, g and U. First, let X := L rn (R). Next, define 



Finally, let 3 be given by 



+oo otherwise 



y(p) = J Jr [dM x ) m/2 )} dx + I'm x2 p( x ) dx if P e W) and d x (p m/2 ) 6 L 2 (R), 
1 +oo otherwise. 

Since any elements p and rj in the proper domain of $ belong to V2 (R), it is clear that = 
g{p,if) = W2(p, rj) implies p — 7/. Further, the lower semi-continuity of J on L m (R) follows 
from Lemma IA.1I in the Appendix [A] Next we show that — for any given c > — the sublevel 
A c := {p € L m (M.)\$(p) < c} is relatively compact in L m (R). To this end, we shall prove below that 
B c := {77 = p m / 2 \p e A c } is relatively compact in L 2 (R); since the map 1 : L 2 (R) L m (R) with 
b{rf) = n 2 l m is continuous, it then follows that A c = l{B c ) is a relatively compact set in L m (M). 
To show relative compactness of B c in L 2 (M), we verify the hypotheses of the Frechet-Kolmogorov 
theorem, see e.g. [T^] Theorem IV. 8. 20]. 

B c is bounded in L 2 (M.). For any given q > 1, the Gagliardo-Nirenberg interpolation inequality 
(in combination with Holder's inequality if m > 2) provides 



f q dx<K^J (d x r)) 2 dx) (^J 77 2 / m dx^ . (69) 

Boundcdness of B c in L 2 (R) follows directly with the choice q = 1; recall that f? 2 /™ = p lies in the 
proper domain of # and thus has integral one. 
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B c is tight under translations. For every r\ £ B c and any h > 0, we have that 



\rj(x + h) — rj(x)\ dx 



d x rj(x + z) dz 



dx < h / (d x r/(x)) dxdz < ch 2 



Thus the integral converges to zero uniformly on B c as h 10. 

Elements of B c are uniformly decaying at infinity. For every rj £ B c and any R > 0, we have 

/ nix) 2 dx < 4 / \x\v(x) 2 dx < ^ ( [ x 2 ri(x) 2/m dx] ( [ n{xf- 2 l m dx) ' . 

J\x\>R R JR H\Jr J \Jr J 

Inside the last expression, the first integral is less than c, and also the second one is controlled in 
terms of c, using inequality (|69p with q = 2 — 1/m > 1. 

In conclusion, B c satisfies the hypotheses of the Frechet-Kolmogorov compactness theorem. It 
follows that J has compact sublevels. 

We turn to verify the two hypotheses in (|68p for the set U := {p Tk \k e N}. The first hypothesis 
is satisfied because of ([60]) and (|6T|). We establish the second hypotheses as a consequence of the 
T-uniform approximate Holder continuity (|67p. For this, pick T > and /i s]0, 1[ arbitrary. Given 
fc e N, define iVfc G N such that (Nk - l)rk < T < Nutu- We distinguish two cases. If < h < Tk, 
then (writing r = Tk and N = Nk for ease of notation) 

r T-h N-l , oo N 1/2 

/ W 2 (M* + &)>/*(*)) dt < r ^ JW,^ 1 ^) < hN^ 2 (J2^2(pi + \pi) 2 



< fc(2riV.F(p°)) 1/2 < (2(T + l)^)) 1/2 / l 
by ([STjl . If, on the other hand, h > Tk, then there is a J € N with h < Jrk < 2h, and so 



L 



T-h N—XJ—X 

W 2 (p T (t + h), p r (t)) d* < T £ ]T ^2{ P n T +j+ \p n T +j ) 



n=0 3=0 

JV-1 , r oo -J l/2x 

<r^ (j 1 / 2 £W 2 (^ +1 ,^) 2 ) <^r(2JrJ-(p )) 1/2 <2(T+l)^(p ) 1 / 2 / l 1 / 2 . 

n=0 ^ L «=0 -I ' 

Theorem 14.91 now provides for every subsequence (tv) C (rfc) the existence of a subsubsequence 
(Tk") Q (rk") such that p T „ converges in measure with respect to t G]0,T[ in L m (M.) to some 
limit p + . By convergence of p Tk (t) to in W2 for every t € [0, T], it follows that p + = p*. By 
the usual arguments, we conclude that the entire sequence (p Tk ) converges to p„ in measure. In 
combination with the r- uniform bound (f5T)f of p T in L m (K), we can invoke Lebesgue's dominated 
convergence theorem to conclude strong convergence of p Tk to in L m (0, T; L m (M.)). □ 

Corollary 4.10. For every T > 0, we have pT /2 £ L 2 (0,T; H\R)) andd x (p™) € L^O, T; L X (R)). 
Proof. Fix T > 0. By estimate p™,^ 2 is uniformly bounded in the reflexive Banach space 



L 2 (0, T; By Alaoglu's theorem, there is a subsequence (tv) C (rk) such that p™^ 2 

converges weakly to some limit £ in that space. Since p" 1 /, 2 converges strongly to p™'^ 2 in 
L 2 (0,T;L 2 (R)) by Proposition (O, it follows that pT /2 = (£ L 2 (0, T; H 1 (R)). 

The second claim is a trivial consequence of the representation d x (p™) = 2p" l ^ 2 d x (pT^ 2 )- □ 

4.5. Derivation of the entropy formulation. In this subsection we show that m* = S _1 [p»] 
satisfies the entropy formulation (|27|) . The following proposition plays the role of Lemma 2.4 in 
[18| . We emphasize that although the entropy inequality contains the same dissipation term as in 
|18j . we derive it from a completely different source. Here, it results naturally from the variational 
construction by minimizing movements. 

Proposition 4.11. Define u* = S -1 [p*] from p„ via scaling, see Then satisfies the entropy 
inequality (f2"T)l for any k > and for any non-negative test function tp € C£°(R+ X R). 
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Proof. In view of Remark 12. 61 it suffices to prove the estimate ([2"T|) for all test functions ip of the 
form ip(t,y) = 6(t)<f)(y), with arbitrary non-negative 9 € C™(R + ) and 4> e C£°(R). Let and 6 
as well as /c > be fixed in the following. 

The proof of the entropy inequality results from another application of Lemma 14.21 For given 
parameters v > and < e < k, define 



>(x) 

where the function S t : R — > R is given by 



(70) 



S e (s) := / Sgn e (r m -fc m )dr. 
Jo 

From its definition, it is obvious that S e (s) is bounded from below by —A: and converges monoton- 
ically (from above) and uniformly in s € R to \s — k\ — k in the limit e J, 0. The first and second 
derivatives are given by 

S' e (s) = Sgn e (s m - k m ), S' e '(s) = 2ms m - 1 6 £ (s m - k m ). 

For later reference, we observe that there is a constant K e such that 

\S e (s)\<s, \S' t {s)\<\, < S'J(s) < Res' 1 forallseM+, (71) 

which can be verified by elementary calculations, using the definitions of 5 e and Sgn e . 

As a preliminary step, we show that there is a K-flow associated to $ which is given as the 
solution operator © e „ to 

d sV = (r) [ Sgn e ( (I) ™ - hT) <f> o T] + ^ (72) 

This will be achieved by application of Lemma 13.41 In the situation at hand, we have 

F(x, V ) = ^(^^oT^T'Or) +1/77 log 77, 

and it is easily seen that (|3"4"|) is satisfied (with c= v > 0). The associated function H reads 

JT(x, = tS e (-^rr) </> o X^X'^) - v log e 



We calculate the entries of the matrix in J37] 



^(^0 = ^(^)/ 1 (x)+^(^)/ 2 (,)+r 1 ^(^)/ 3 (,), 
^(«.0 = 5.(^)/^)+r^(^)M*)+r^'(^)/.(x). 



Here /i to fe are continuous functions of compact support in R, explicitly expressible in terms of 
X and its derivatives. Using the properties (|7lj) . it is easily seen that there is a constant M e such 
that 

holds for all £ G R+ and uniformly in x £ R. Thus, for a suitable choice of k = K e ,i/, the matrix 
in (|37|) is positive for every £ £ R + and a; € R, and so H is jointly convex. By Lemma \3. 41 the 
solution operator 6 £i „ for (|T2[) is a K £: „-flow for \& ejJX . 

In order to apply Corollary 14.31 we need to calculate the derivative of T along solutions 77 
to ([72")) . For simplification, introduce the rescaling v(t) = S _1 [77(t)], see ([SJ) . Further, recall the 
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definitions of 1 and of a — 1' o 1 1 , and properties (|19l) and (j^Uj) . Using the functional defined 
in (|65[) . we have 

as Jm Jr 

<- f {a«}o^[{(aoX- 1 )-Hr- 1 }ol] [{Sgn £ (« m -P)^}oS] da: - i/Jfcfa) 

f ^J^^V 1 " 1 ] [Sgn e (v m -k m )<j>] dy-vR^r,) 
J R Lm— 1 Jy L J f 

/ + &«™] [Sgn e (u m - F")# y dy - i/jfefo) 

Jr 

= - jf [(« m )„ + 6(« m -A; m )] Sgn e ( W m -fc"%,dy- / Sgn e (tT - y dy 

[{v m ) y + b(v m - k m )] (v m ) v Sgn' e (v m - k m )<f>dy - i/Jfefa) 

= - / <^,S e («) dy+ f (b(v m - k m )$ v - k m b y 4>) Sgn e (v m - k m ) dy 

[P e (v) y ] <f>dy - I bQ e (v) y (j)dy - v&h(t)). 
Jm 

In the last step, we have implicitly defined the smooth functions P e ,Q e : R+ — > R such that 

P' t {s) 2 = Sgn' e (s m - k m ) and Q' £ («) = (s m - k m ) Sgn' e (s m - k m ). 
Accordingly, we define, still with v = S _1 [r/], 

&cAv) = - I <t>yySM dy + [ (b(v m - k m )<^ y - k m b y 4>) Sgn £ (v m - k m ) dy 
Js. Jm 

+ f [P^vU^dy- [(b^yQ^dy + vRHiv)- 



Considered as a functional of v, the right-hand side is lower semi-continuous with respect to strong 
convergence of v in i Tn (R). Indeed, all of the integral expressions are even continuous in L m (R), 
except for the one involving P e , for which lower semi-continuity can be concluded by means of 
Lemma [A. 11 Since convergence of r\ in L m (M.) is equivalent to convergence of v = S -1 ^] ini m (R), 
the functional & e .„ is lower-semicontinuous with respect to rj. Now since 6f v r\§ converges to ?/o m 
L m (R) for every t\q <E L m (M.), we conclude that K eM satisfies the condition ([57]) for the application 
of Corollary [OJ 

We shall now derive a refined version of estimate (|59|) . To this end, recall that 8 <G Cj?°(R+) is 
a temporal test function. Multiply (|58|l by 6 (tit) and sum over n to find 

r ± ^Ap n r) e{nr) ~ 6 l in + 1)T) > r ± ^(PrMnr) + ^(p°)- 

n—l n—1 

Thanks to the strong convergence of p T to in L m (R) and the lower semi- continuity of A e ^, we 
can pass to the time-continuous limit and find 

* e ,„O>.(t))0'(t)dt> / 6{t)& f _. v { P *{t))&t. 



Taking into account estimate (|63p on H(p*(t)), we can now pass to the inviscid limit v 1 and 
find, recalling f(t,y) = 9(t)<f)(y), that 

f T [ S t (u* - k m )d t (pdydt > f T [ ([(<% + 6« l - k m )]ip y - k m b yV ) Sgn £ «" - k m )dydt 
Jo Js. Jo Jm v ' 

+ f [ [P t (u*) y ] 2 <pdydt- [ [ (b(p)yQ e (u*)dydt. 
Jo Jm. Jo Jm 
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Notice the integration by parts in the first term, which is admissible since d x (p™) € L 1 (0, T; L 1 (R)) 
by Corollary 14.101 In the final step, we pass to the limit e j 0. By the uniform convergence of 
S e (s) to \s — k\ — k, by the uniform convergence of Q t to zero, and since 

[PM y Y =Sgnl £ (u?-k m )[(uT)y] 2 

by definition of P € , we finally obtain ([271) . □ 



4.6. L°° bounds. Definition 12.41 of entropy solutions requires u* G L 1 n L°°(]0,T[xK) for every 
T > 0. The Z^-bound is obvious from our construction. Below, we prove the Z/°°-bound. 

Proposition 4.12. Assume that p" G 7> 2 ( R ) n L°°(R), and let k > be such that p°(x) < 
ka{x)" 1 / {m ~ 1 "> for almost all x G R. T/ien p*(T,x) < ka(x)- 1 / { - m -^ /or all T > and almost 
every i£l. 



Proof. Once again, the proof is obtained in application of Lemma 14.21 As auxiliary functional, we 
choose 

^e,v(v) = / PosPar e {a(x)~^ r]{x) - fc)a(x)™^ 5 ' dx + vW{r\) 

with f > and e £]0,k[. We verify that \l/ £ .„ satisfies the assumptions required in Lemma 13.41 
With the short hand notation A = a x ^ m ~ 1 ', we have 

'A(x) 

Recall that PosPar^ = Stp e and PosPar" = S ei and observe that there is some e such that 

< PosPar e (s - k) < s, < Stp e (s - ft) < 1, < 5 e (s - k) < K^s' 1 (73) 
for all s G R+. We obtain for the second derivatives of H: 
A{x) fA(x) 



H(x,£) = £ PosPar, (J^p- - fc)A(x)- 1 - v\og^. 



H ii {x,i) = ^5 e (^--k)+^-\ 



H xx (x,£) = /i(x)^PosPar £ - k) + f 2 (x) St Pe - k) + h(x) * 



5, -^--fc 



Hxd x >0 = /4(a;)PosPar e ^— fcj + — — Stp e ^— fcj + -J^.^— 1 ft), 

where the functions f\ to fa are explicitly expressible in terms of A, A' and A" , and are uniformly 
bounded in x G R. Taking into account (|T3[) and that A > a 1 /^™ -1 ), we arrive at the uniform 
bounds 

with a suitable constant M e . Consequently, by application of Lemma f3.4[ the solution operator 
(5<ii to 

d s r) = B x (r]B x [Stp e (a ^i?/ - ft)]) + 

defines a K-flow for ^e^. The regularizing effect of the viscous term is strong enough to justify 
the following calculations: 

-^.F(6^ ) = j a^- 1 ( v [St Pe (a^ V - k )] x ) x dx + v J arf 1 " V* dx 

<-j ^[(a^)™- 1 ] JStp £ (a^ - k)] x dx - ^ dx. 

Here we have used that the second integral with v can be estimated as in (|S4"|). neglecting a 
positive term. The product under the first integral is always non-negative since both the functions 
s i-> s m_1 and s H> Stp £ (s — k) arc differentiable and increasing. Thus, condition ([57| is satisfied 
with 



Jtfo) := 1± / n m dx 
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which is obviously continuous with respect to strong convergence in L m (IR). Estimate (|59l) in 
combination with (|50[) yields 

N r 

<*Ap a )+( — NT + 2T(-K))F(p°). 

\ ma 1 

For fixed positive parameters v and e, the modulus k of convexity is a r-independcnt constant. 
We can thus pass to the limit Tk i in (|59|) and obtain, using the lower semi- continuity of ^S e . v in 
W 2 , that 

Mp*(t))<Mp°) + —tf( p °) 

ma 

for every T > 0. Using further that H(/9*(T)) is a finite quantity, sec (|63|) . we can pass to the 
limit v 4- and obtain 

/ PosPar e (a(x)^ rT p*(T, x) - fc)a(x)^=r dx < / PosPar e (a(x)'^p (x) - fc)a(a;)^ T dx. 

By the properties of a, and since p°,p*(t) G i m (M), we can further pass to the limit e l 0, which 
yields 

/ [p»(r,x) - fca(x)^] dx < / [p°(x)-fca(x)^r] dx. 
it Jr 

Since the integral on the right-hand side is zero by hypothesis, so the is the integral on the left-hand 
side. This proves the claim. □ 

Corollary 4.13. Provided that p° e L°°(R), it follows that p™ e L 2 (0, T; _ff 1 (R)). 

Proof. By Proposition 14. 121 

K := sup esssup,o*(i,x) < oo, 

teR + xeR 

and so 

" [dMt) m )} 2 dx = [ [2 P ^t) m ' 2 dMt) m ' 2 )] 2 dx<AK m [ [dMt)™/ 2 )} 2 dx 

JR JR 

for alH > 0. The claim now follows from Corollary 14. 101 □ 
4.7. Continuity at t = 0. It remains to verify that it* attains the initial condition u°. 
Proposition 4.14. u*(*) -> u" in L 1 (R) as < | 0. 

Proof. Since the rescaling S from ([5]) is a homeomorphism on L 1 (K), it suffices to show that 
p*(t) — > p° in L 1 (K) as £ | 0. By the lower semi- continuity of J-", we can pass to the time- 
continuous limit t ^ in (|50[) and obtain 

•F(p.(*)) < liminf HPr{t)) < Hp°) 

tJ.0 

for every t > 0, which implies that 

limsup ?{p*(t)) < Tip )- (74) 

On the other hand, since the limiting curve is continuous in W 2 , once again the lower semi- 
continuity of T in W2 yields 

T{p°) < Hminf T{p*{t)). 

In combination with (|74p . we have that 

M%(t))=J(;). (75) 
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By definition, J-(p) is the mth power of the L m -norm of p with respect to the non-uniform 
background measure m~ 1 a(x) dx. Since m > 1, and since the weight function a satisfies the 
bounds (al)&(a2), we can conclude by standard arguments, see e.g. [21] Theorem 2.11], that 
weak convergence of p*(t) and convergence (|75)l together imply strong convergence of p*(t) to 
p° in L m (M.). To obtain convergence in L 1 (IR), we apply the generalized Holder inequality with 
exponents 2m/(m — 1), 2 and 2m: 



f \ P *(t)-p°\dx< [ (l + \x\ 2 )-?[(l + \x\ 2 )\(p4t)+p )^[\p*(t)-p°\^dx 
Jm. Jr. 

(l + \x\ 2 )'^ dx 



(l + \xn(p*(t)+p»)dx 



\\p*(t)-p°\\ 



The first integral on the right hand side is clearly finite, and the second integral remains uniformly 
bounded as t ], 0, since continuity in W2 implies continuity of the second moment. The last term 
vanishes for t J, because of the strong convergence of p*(t) in L m (M.). □ 



4.8. Proof of Theorem 12. Tl At this point, we have proven that it* meets all the requirements 
for being an entropy solution as stated in Definition ^. 41 we have a, Ei'fl L°°(K) by Proposition 
EE] and it™ G L 2 (0, T; _ff 1 (R)) by Corollary HT3] we further have continuity of at t = 
by Lemma 14.141 and the validity of the entropy inequality has been verified in Proposition 14. 1 fl 
Finally, convergence of the time-discrete approximation scheme in Wasserstein and in L m have 
been shown in Lemma 14.71 and Proposition 14.81 respectively. 



5. Uniqueness of entropy solutions 



In this section we prove Thcorcm l2.8[ using the doubling of the variables device. Since we follow 
almost literally the proof of [TTJ Theorem 1.1], we restrict ourselves to the key steps and refer the 
interested reader to the original article 11 for more details. 

For a fixed T > wc shall use the notation LI T = Kx]0,T[. Let <p e C^°(Ut x Ht) be a non- 
negative test function, and assume that u and v are entropy solutions in the sense of Definition 
2.41 For brevity, we write v = v (x, t), u = u{y, s) and (p — cp(x, t, y, s). From Proposition 14. 1 ll wc 
obtain 



(|« - u\<pt + Sgn(« - u)[(v m - u m )b(x) - (v m ) x ]<p x 

n T JJn T 

- Sgn(u - u)b'(x)u m ip) dx dt dy ds (76) 

<-limsup/7 ff (v m ) 2 x Sgn' e (v m - u m )ipdxdtdyds. 
eio JJiit JJn T 

Recall that Sgn e is a smooth uniformly convergent approximation of the sign function, which is 
obtained by mollification with 5 e . Since (v m ) x G L 2 (TLt), the following integration by parts is 
justified: 

- [f Sgn e (v m - u m )(v m )*Vy dyds = [[ (Sga e (v m - u m )) y (v m ) x <pdyds. (77) 
JJn T JJn T 

We integrate ((77)) w.r.t. (x, t) G ±1 T and send e \. 0. By the dominated convergence theorem, that 
yields 

Sgn(i> — u)(v m ) x (p y dy ds dx dt 

(78) 



Mm // (u m ) y {v m ) x Sgn[(v m - u m )ip dy ds dx dt. 
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Adding (JTSJ) and ([78]) we get 

(\v - u\<p t + Sgn(v - u)[(v m - u m )b(x)ip x - (v m ) x (ip x + tp y )] 



— Sgn(u — u)b'(x)u m ip) dxdtdyds 



(79) 



< — lim 



imsup/7 // [(v m )l-(u m )y(v n %]Sg<(v 
e^O JJn T JJn T 



m „ , m 



u m )ipdx dt dy ds. 



The terms in (|T9[) containing b can be rewritten as follows: 

Sgn(w - u)(v m - u m )b{x)<p x - Sgn(v - u)b'(x)u m <p 
= Sgn(v - u){v m b{x) - u m b(y))^ x + Sgn(v - u)(u m (b(y) - b(x))<p) a 



(80) 



Now, we repeat the previous steps with a simultaneous interchange of the roles of u and v and the 
roles of (y, s) and (x,t). Summation of (|79j) with its respective counter part yields 



< 



(|t> - u\(tpt + (p.) + Sgn(« - u) [v m b{x) - u m b{y)} (<p x + ip y ) 
n T JJn T 

+ \l! m - U m \(cp xX + 2^ + iPyy) 

+ Sgn(v - u)[{u m (b(y) - b(x))<p) x - (v m (b(x) - b(y))<p) v ] ) Ax dt dy ds 
-hmsup if [[ \{v m ) x - {u m ) y \ 2 Sgn' e (w m - u m )(p dxdtdyds < 0, 



(81) 



where we have used the identity 

Sgn(v - u){{v m ) x - (u m ) y ) = [\v m - u m \ x + \v m - u m \ y ] 

and integrated by parts. We emphasize that for the estimation on the right-hand side of (|8ip . the 
dissipation term from the definition (j2"7| has been essential. 
As usual, tp = ip{x, t, y, s) is chosen in product form, 

tp[X,t,y,s) = i>{—^—, — Y~ )Ucr{— — )5a{—^- j, 

where < tp £ C%°(ILt) is a test function, 8 a = o~ 1 5i(t/<j) for a > and Si is the mollificr 
defined in (0, and 



u a (x) 

Accordingly, we introduce a new set of variables (x,i,z,r) with 



x + y - t + s 
x = — - — , t 



x — y t — s 

z = — - — , r 



2 ' " 2 ' 2 
for which we have (by the usual abuse of notation) 

dt=d t + d s , d x = d x + d y , d xx = d xx + d. 



2 ' 



, 28 

yy i &u xy . 



With the understanding that u = u(y,s) and v = v(x,t) while ip = ip(x,i), u a = u a (z) and 
8 a = 8 a (r) 1 the inequality (|81l) can be written as 



0>Ja 



n T JJn 



[\v- u\dpp + Sgn(« - u){v m b{x) - u m b{y))d x tp + \v m ~ u m \d xx tp 

f Sgn(v - u)[((u m (b(y) - b(x))) x - (v m (b(x) b(y))) v )ip 
+ u m (b(y) - b{x))d x il) - v m {b(x) - b{y))d y if\ }oj a S a 



(82) 



(b(x) -b(y))\v m - u m \ip8 a d,u) 



dx dt dy ds. 
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By assumption, b £ W 1,oa (M.). Moreover, since u and v belong to L°°(Ht), and since S[(x) < 
for x > 0, we can find a constant K (depending on b) such that 

(6(x) - b(y))\v m - u m \d z u a = (b(x) - b(y))\v m - u m \zo-- 3 6[(^) 

(83) 

z 

— ' 

T — 



X 

< K\v - u\—^a~X\z\<2a- 



Now we perform the limit a \. 0, which concentrates the support of tp on the diagonals x = y and 
t = s. By Lebesgue differentiation theorem we then obtain 



> lim J a > - 



\\v- u\%Ih + b{x)\v(x, t) m - u(x, t) m \^ x + \v(x, t) m - u(x, t) m \^ xx 
+ b'{x)\v(x,t) m -u(x,t) m \4>} + K\v(x,t) -u(x,t)\ip 



(84) 



dx dt. 



Regrouping terms and using again that b £ W 1 '°°(M.), we arrive at the key estimate 

(\u-v\4) t + b\v m ~u m \4j x + \v m -u m \ip xx )Axdt< C [[ \v-u\ipdxdt, (85) 
n T JJn T 

for some C > depending on K and on b. 

At this point, we make the classical special choice for the test function ip. For any given 
< ti < £2 < T and r > 0, let a £ C£°(]0, T[) be the cr-mollification of the characteristic function 
of the interval [tijta], and let 4» r £ C^°(R) be such that c/> r (x) = 1 for \x\ < r and 4> r (x) = for 
\x\ > r+ 1. Setting ip(x,t) = <j> r (x)8 cr (t) yields 

lim ff (b\v m - u m \4> x + \v m - u m \ip xx ) dx dt 



r— >oc 



Hi 



r~ >oo 



< C lim // (v + u) dx dt = 



by the dominated convergence theorem (with respect to t £ [0, T]), because u and v arc t- uniformly 
bounded in L 1 (M). Therefore, by sending r — > 00 in (f85|) . we get 

-[ [ \v-u\e' a (t)dxdt < C [ ( \v-u\6 a (t)dxdt. (86) 
Jo Jr Jo Jr 

Finally, passing to a J, 0, we obtain 

IK*a) - v{t 2 )\\ Ll < \\uih) - v(Jt!)\\ L i +C [' \\u{t) - v{t)\\v dr, (87) 



hence we can use the Gronwall inequality in (|87[) (integral form for measurable functions, cf. |13j ) 
to obtain 

\\u(t 2 )-v(t 2 )\\ Ll < IKtO-vCiiJUiiCl + C^a-tOe ^-* 1 )). (88) 

Since u and v are right continuous at t = 0, we can perform the limit t\ \. in (|88[) and obtain 
stability of the entropy solutions u, u in the i 1 (R) norm. In particular, if v(0) = u(0), then 
v(t) = u(t) for all t > 0, which shows uniqueness. 

Appendix A. A lemma on lower semi-continuity 

Lemma A.l. Let H £ C 1 (]R+), and let <f> £ C°(R) be a bounded non-negative function. Define 
the functional * : V 2 (R) -4 [0, 00] by 

*(f/) = [ 4>(x)[d x H( V (x))] 2 dx 



whenever the integral is well-defined, and +00 otherwise. Then is sequentially lower semi- 
continuous in L P (JS>), for arbitrary p > 1. 
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Proof. Let (r/ n )„ e N be a sequence that converges to T]q in L P (R), with 'J' := sup„ v I'(ry„) < oo. 
Without loss of generality, we may even assume that r] n converges to rj almost everywhere on K. 
For every e > 0, introduce the functional "J,, : ^(R) ~> [0, oo] with 

\& e (rj)= I <j)[d x H(ri)] 2 Ax where Q c := {x e [-e"\ e~ l ]\(f>{x) > e}, 

with the understanding that ^(j?) = +oo unless H(rj) £ /f 1 (f2 c ). Since <f> is non-negative, 

e I [d x H(r] n )} 2 dx < < *(r?„) < ¥. 

Thus, the functions H{r) n ) are n-uniformly bounded in H 1 (f2 e ). By Alaoglu's theorem, every subse- 
quence (H(r/ n r)) n i£f{ contains a subsubsequence (H (7?n"))n"eN that converges weakly in _ff 1 (f2 e ) to 
some limit /i. Further, since fi e is compact by definition, Rcllich's Lemma applies, and H{rj n ii) con- 
verges to h strongly in L 2 (f2 e ). We have assumed that r\ n converges pointwise almost everywhere 
to 770, hence H(r] n ») converges to H(r] ) almost everywhere, and we conclude h = H(r]o) € lf 1 (0 € ) 
— independently of the chosen subsequence. From here, it follows that 

d x H{r ln )^d x H( m ) in£ 2 (fi £ ). (89) 

By elementary calculations, one verifies that 

*e(Vn) > *M + 2 f 4>{p x H{r, n ) - d x H( m ))d x H( m ) dx. 

Using (|59^) . the limit n — > 00 provides 

^e{Vo) < liminf * e (ry„) < 

n—too 

To conclude the proof, observe that 4' e (?7o) - > ^(vo) m the limit e 1 by the monotone convergence 
theorem. □ 
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